arXiv:1502.02684vl [quant-ph] 9 Feb 2015 


Universal two-qubit interactions, measurement and cooling for quantum simulation 
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By coupling pairs of superconducting qubits through a small Josephson junction with a time- 
dependent flux bias, we show that arbitrary interactions involving any combination of Pauli matrices 
can be generated with a small number of drive tones applied through the flux bias of the coupling 
junction. We then demonstrate that similar (though not fully universal) results can be achieved 
in capacitively coupled qubits by exploiting the higher energy states of the devices through multi¬ 
photon drive signals applied to the qubits’ flux degrees of freedom. By using this mechanism to 
couple a qubit to a detuned resonator, the qubit’s rotating frame state can be non-destructively mea¬ 
sured along any direction on the Bloch sphere. Finally, we describe how the frequency-converting 
nature of the couplings can be used to engineer a mechanism analogous to dynamic nuclear polar¬ 
ization in NMR systems, capable of cooling an array of qubits well below the ambient temperature, 
and outline how higher order interactions, such as local 3-body terms, can be engineered through 
the same couplings. Our results demonstrate that a programmable quantum simulator for large 
classes of interacting spin models could be engineered with the same physical hardware. 

PACS numbers: 


I. INTRODUCTION 

Quantum computing promises revolutionary improve¬ 
ments over traditional classical algorithms to solve hard 
problems, provided that a functional, noise-tolerant 
quantum computer can be built. A number of paradigms 
for engineering a quantum computer have been consid¬ 
ered in recent years, with the most studied being the 
“gate model,” where the Hamiltonian of a system of N 
quantum bits (two-level systems) is assumed to be triv¬ 
ial and all unitary evolution of the system wavefunction 
occurs through a digital series of finite-length pulses to 
enact specific unitary operators. A key fact which makes 
such constructions feasible is that the sequential combi¬ 
nation of single-qubit rotations with a two-qubit inter¬ 
action suffices to generate all possible two-qubit unitary 
transformations [T]. This allows arbitrary quantum op¬ 
erations to be performed through a small set of distinct 
physical terms, greatly simplifying the underlying system 
architecture. 

However, in many cases it is desirable (or unavoidable) 
for the system to have a nontrivial physical Hamiltonian 
H that contributes continuously to its time evolution, 
and since single-qubit rotations generically do not com¬ 
mute with H, performing arbitrary deformations of the 
system wavefunction through a small handful of discrete 
gates becomes difficult or impossible. Further, in the adi¬ 
abatic model of quantum computing [5] , only slow varia¬ 
tions of the system Hamiltonian are allowed and thus all 
unitary time evolution must occur from the application of 
a continuous, physical H (t). Quantum computing with 
physical Hamiltonians has a number of potential advan¬ 
tages over the pure gate model, including simplicity of the 
algorithms |3] , resilience to low frequency noise and small 
variations in the system parameters (often through topo¬ 
logical protection HIS]), and even passive correction of 


quantum errors arising from unwanted interactions with 
the environment [SHni- Further, such Hamiltonians could 
be tuned to mimic the structure of interesting spin mod¬ 
els in condensed matter physics, directly simulating ex¬ 
otic condensed matter systems in a similar fashion to cold 
atoms in optical lattices m- A simple method for con¬ 
tinuously generating arbitrary multi-qubit Hamiltonians 
is thus extremely desirable. 


In this article, we demonstrate that in circuit QED 
systems of superconducting qubits, a Josephson junction 
coupling with a time-dependent flux bias serves as a uni¬ 
versal interaction resource, and can generate arbitrary 
qubit Hamiltonians of the form CabO'iCr^, where 

the nine Cab coefficients generate all possible two-body 
interactions. We then show that by capacitively cou¬ 
pling a qubit to a resonator and applying appropriate 
multi-photon signals through the qubit’s flux line, we 
can non-destructively measure the rotating frame state 
of the qubit along an arbitrary direction on the Bloch 
sphere. This technique does not require that the many- 
body Hamiltonian is turned off during the measurement, 
allowing aspects of the system to be probed without dis¬ 
rupting the many-body state. We also propose a sim¬ 
ple mechanism based on these couplings which can cool 
a grid of qubits far below the ambient temperature of 
their environment. Appropriately tuned, this scheme 
could be used to aid the operation of an adiabatic quan¬ 
tum computer. Finally, we demonstrate that frequency¬ 
converting couplings (which promote a tunneling photon 
up or down in energy as it hops between qubits) can 
be used to effectively eliminate the qubit nonlinearities, 
generating higher-order interactions than simple 2-body 
terms. 
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FIG. 1: Basic coupling types considered in this work. In 
figure (a), we demonstrate a flux biased Josephson junction 
coupling between two transmon qubits (blue dashed regions), 
discussed in Sec. [TI] each of which has its primary Josephson 
junction split into a SQUID loop for optional flux tuning. By 
tuning the flux signal F {t) threaded through the loop made 
by the smaller junction and the bridged ground, all possible 
two-body interactions between the qubits can be dynamically 
generated. Similar results can be derived for flux and flux- 
onium qubits in the appropriate parameter regimes. This 
cou pling can be used to implement the cooling protocol of 
Sec. IV In (b) we depict a transmon qubit coupled to a res¬ 


onator through a simple capacitive coupling, as discussed in 
Sec. |III| The grounds are bridged in this figure, though unlike 
case (a) this is not required to obtain the desired behavior. 
A flux signal is applied through the transmon’s SQUID loop 
to generate couplings which allow for the transmon’s rotating 
frame state to be measured along any direction on the Bloch 
sphere. This setup can also be used to couple qubits, though 
without as much flexibility as the small junction in (a). 


II. UNIVERSAL TWO-QUBIT HAMILTONIANS 

We consider a pair of charge insensitive superconduct¬ 
ing qubits, such as transmon, flux or fluxonium devices 
[min]. These qubits have one or more quantum phase 
degrees of freedom (j), and when operated at a flux sym¬ 
metry point and restricted to the lowest two levels, the 
phase operators take the generic form 


they share the same physical ground and are otherwise 
coupled by a Josephson junction of strength aEj (where 
Ej is the internal Josephson energy of the qubits), as 
first studied in [Ulli], with a time dependent flux bias 
F (t). The grounds of the two circuits are bridged so that 
charge may flow through the devices and coupling junc¬ 
tion in a closed loop, as in FIG. 1. We assume that the 
flux is solely threaded through the coupling region and 
does not leak into any internal loops in the qubit devices, 
and Anally assume that the inductances created by bridg¬ 
ing the grounds are negligible. Appealing to the circuit 
quantization equations for superconducting devices m, 
the interaction term takes the form 

Hint = -aEj cos {(1)1 - (j)2- F (t)). (2) 

We assume that a is small so that in the absence of 
any applied drive tones the qubits do not interact and 
the distortion in the qubits’ energy levels is minimal M- 
The bare excitation energies of the qubits are wi and ui 2 ■ 
Anticipating that we will drive the system through F (t ), 
we will work in a resonantly rotating frame through the 
transformation Ifk) —)• exp [i + ‘^^2) A 

in turn adds time dependence to the qubit operators: 

^ af ^ ( 3 ) 

We now choose: 

F {t) = 7r/2 - f ^^-2 cos [(wi +uj2)t + xi+ X2J 
+fx°J cos [(wi -UJ2)t + Xl- X2]) ( 4 ) 

+2y2(/..cos^^^i^ 

Here, the x ^nd '0 terms are constant phase shifts, the 
coefficient /3 controls the relative magnitude of the 2- 
photon and frequency-converting terms, and we assume 
all four of fzz,fxy,fxz and fzx are small compared to 
unity. Expanding the interaction term, we have 


+ fzx cos 


UJ2t J- 'lj)2 


sin (j) ^ s {a^ + a ) , cos </)—>■ cq J- ctr^. (I) 

The charge operator Q becomes and z is the natu¬ 
ral rest frame axis of the qubit. For devices such as flux 
qubits with two or more degrees of freedom, the phase 
operators can still nearly always be written in this form, 
with additional coefficients and sign flips if necessary. We 
now couple two qubits together, wiring them such that 


Hint = —ctEj [(cos 01 cos 02 + sin 0i sin 02) cos F (t) 

+ (cos 02 sin 01 — cos 0i sin 02) sin F (f)]. (5) 

We will now expand the terms in Hint to lowest nontriv¬ 
ial order in the /’s, keeping only the terms which are 
independent of time and neglecting everything which is 
rapidly oscillating. The result for is: 
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Signal 

Rotating frame coupling 

Role 

fzz (Constant bias) 
cos [(tJi -bJ2)t + Xl-X2\ 
cos [(tJi - 1 - 0^2) t + Xi + X2] 
Uz cos 

(co -|- CUi) (co J- C(t|) 
g-i(xi-X2)^+^- gi(xi-X2)^-g,+ 

g-*(xi+X2)g,-+g-+ _|_ g»(xi-l-X2)g--g,-- 

(g-iV>l^+ + (g^ + g^l) 

2:2; interaction 

Complex photon hopping/frequency conversion [ 16 ] 
2 -photon coherent pump 
State-sensitive pump (continuous CNOT) 


TABLE I: Fundamental interactions generated from driving a flux-biased Josephson junction with the signal Q. As shown in 
the text and equation the nine free parameters of the signal allow for all two-body Hamiltonians (J2 to be 

generated from the same fixed device layout, FIG. la. Single qubit « terms can be eliminated by detuning the drive signal(s) 
or applying flux biases to the qubits’ internal flux lines; x and y terms can be eliminated by driving the qubits individually. 


Hint = -aEj f^^{co + cai){co + ca2) + fi°J {cos XiUi + sin xicrf) {cos X2cr2 + sin X2cr2) 

+fL (cos V'lCTi -b sintpiaf) (cq + ccr^) + (cosV'2CT2 + sin?/)2crf) (co + caf)] + O (f) 

= CabCTj(T2 + (single qubit terms). 


( 6 ) 


The single qubit terms can be cancelled out by com¬ 
binations of detunings and single photon drive fields, if 
desired. Counting degrees of freedom (including the rela¬ 
tive magnitude of the 2-photon and frequency-converting 
parts of the fxy signal, which we set equal to 1 to make 
the above expression more compact), we see that we have 
nine independently adjustable parameters, sufficient to 
independently tune all nine of the Cab- We note that in 
the more general case where the qubits have internal flux 
biases which are not 0 or tt (and time reversal symme¬ 
try is thus broken), instead of using the trigonometric 
operator identifications Q it can be more convenient 
to expand the coupling in the qubit phases by letting 
(j) = 4>o + a {< 7 '^ + cr~), where (jjQ is the minimum of the 
biased Josephson potential in each qubit. One will arrive 
at the same result in either case. 

We have thus demonstrated that a flux biased Joseph¬ 
son junction is a truly universal interaction resource 
for quantum Hamiltonian engineering in superconduct¬ 
ing devices. A wide array of interesting spin models 
can be implemented through these couplings- for exam¬ 
ple, Kitaev’s honeycomb m model can easily be en¬ 
gineered in both the abelian and non-abelian parame¬ 
ter regimes, though the question of how to prepare and 
maintain the ground state in either case is subtle and 
difficult. We note also that neglecting rapidly oscillating 
terms requires that all matrix elements are small com¬ 
pared to tci/2, UJ 2 / 2 , Iwi — W 2 I, uji + UJ 2 and so forth. 
Consequently, the qubit detunings and nonlinearities set 
the ultimate limit on the strength of interactions which 
can be generated through this method. A second limit 
is that continuous application of microwave helds will 
slowly heat the device by creating quasiparticle excita¬ 
tions; however, as the Cooper pair breaking energy in 
commonly used superconductors such as aluminum or ti¬ 


tanium nitride is more than an order of magnitude larger 
than typical drive frequencies in this setup, quasiparticle 
generation will be suppressed. 

As a final note, as first pointed out in [181 IH]) the 
frequency-converting couplings [16j (which conserve to¬ 
tal 2 spin and tunnel a photon from a qubit to one of 
its neighbors) can be used to engineer an artificial gauge 
field, where the electrically neutral photon excitations of 
the qubits behave as if they were charged particles mov¬ 
ing in a magnetic field. Couplings produced by drive sig¬ 
nals with phase delays (the xi X2 phases in Eq. 
become complex in the rotating frame, and if a closed 
loop is formed by couplings where the phases x ho not 
sum to an integer multiple of 27r, then an artificial gauge 
field exists and the loop is pierced by an effective mag¬ 
netic flux. This can be used to engineer single photon 
circulators [20]? and in large lattices, fractional quantum 
Hall states of light m- 


III. ALTERNATIVE FORMULATION FOR 
CAPACITIVELY COUPLED QUBITS: 

ARBITRARY NON-DESTRUCTIVE 
MEASUREMENTS IN THE ROTATING FRAME 

The ability to engineer arbitrary interactions in the ro¬ 
tating frame is of limited usefulness unless it is coupled 
with an equivalent capacity for making arbitrary mea¬ 
surements of the superconducting qubits. In few-qubit 
experiments (or in purely gate-based quantum comput¬ 
ing), arbitrary measurements can always be performed 
by simply turning off all the couplings at a chosen point 
in time and then using single-qubit rotations to rotate 
the qubits onto a natural measurement axis (conven¬ 
tionally, a^) before measuring them through standard 
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means. However, in larger systems, the requirement of 
turning off the qubit Hamiltonian is a severe constraint, 
and is particularly destructive to continuous error cor¬ 
rection schemes (which are generally sensitive to the en¬ 
ergy levels of the many-body Hamiltonian in such sys¬ 
tems). We therefore present a simple mechanism for 
non-destructively measuring the rotating frame state of 
single qubits, using a capacitive coupling to a resonator. 
Of course, one can achieve the same thing through a flux 
biased Josephson junction; however, capacitive couplings 
are substantially simpler to engineer (and do not require 
the grounds to be bridged). We will do so by exploit¬ 
ing the higher excited states of the transmon qubit; see 
for example [ 22 U 24 ] for previous techniques which exploit 
this structure. 

As before, we consider a superconducting qubit T, 
which we will treat as a three level system with non¬ 
linearity — 6 . For transmon qubits, S is positive and 
roughly equal to the charging energy Ec', for flux and 
fluxonium qubits, S is negative and substantially larger. 
The qubit is capacitively coupled to a resonator R, and 
our goal is to engineer the rotating frame Hamiltonian 

H' = Ah • aT + aiij ■ By measuring the voltage sig¬ 
nal at the resonator, we can thus measure the dressed 
state of the qubit. We let the strength of the capaci¬ 
tive coupling between the qubit and resonator be g, and 
let the detuning a;_R — wt = A be large. The system’s 
Hamiltonian is 


H = —-a^ayaTOT ( 7 ) 


+5 {alpaR ■ 


Before applying any drive fields, we want to diagonalize 


this Hamiltonian, for up to two photons in the transmon 
and {n, n -|- 1, n — 1} in the resonator, to lowest order in 
^/A. In the occupation basis \nT^ nji), to lowest order in 
g/A we obtain 


| 0,'^ ) = |0,n) - £ |l,n - 1) v^; (8) 

\hJl) = ~ |2,u- 1) y/n 

+ -^ |0, n -I- 1) -s/n -I- 1; e nuj r + cut■ 

Here, <72 is the bose enhancement of the matrix element 
to add a second photon to the qubit (<72 —>■ V 2 in the limit 
of zero nonlinearity). We ignore | 2 ,n) as we are going 
to perturbatively eliminate the second level to reduce to 
the spin basis. We now apply the following combination 
of drive fields through the transmon’s flux line (quantum 
phase variable (p): 

Hd = 2 cos(j){ficos[{u}R + ojT)t + x] +f2Cos[At-x]) 
-1-2/3 sin (p cos [ujRt]. ( 9 ) 

We now make the operator identifications: 

cosp —)■ Co-I-cOyOT + C2 (|n-I- 2 , m) (n, m|-I-H.c.) ( 10 ) 

sin(() —)• Si (| 1 , m) ( 0 , to| -I-H.c.)-I-S2 (| 2 , to) ( 1 , to| -|-H.c.) 
The explicit flux signal required to generate Hd de¬ 
pends on the details of the qubit. For example, in 
a split transmon | 25 j with parallel Josephson junctions 
Eji and Ej2, and no static flux bias, starting from 
Hd = -Eji cos {(p- (t) /2) - Ej2 cos ((^i + $ (t) /2) we 

get (all three of fci, fe, fca are assumed to be small): 


Hd 


ki cos 


{ujR + ujt) t + x 


k2 cos 


{ujR -UJT)t-X 


^3 cos UJrI 


2 2 

— {Eji + Ej2) cos P + {Eji — Ej 2) ks sin p cos ojRt + (off resonant) 
r 7.2 1 jL2 I jL2 7,2 7.2 


( 11 ) 


{Eji +Ej2)cosp 


- TT COS ((wfl + Wt) t + x) - — cos {{ujR - Wt) t - x) 


0 (fc 4 


One can then read off the coefficients /i, /2, /a in Q from 
( 11 ). We also obtain a term proportional of the form 
^3 cos p cos 2 u!Rt; however, this term will only participate 
at higher orders in perturbation theory and will have a 
negligible impact. 


note that 


(l,n — l| cosp |0,n) 
(l,n -I- l| cosp | 0 ,n) 


We now move to the resonantly rotating frame where 
the qubit rotates at lot and the resonator at lor. We first 


(0, n -|- l| sin(() |0, n) 
(l,n -|- l| sin/) |l,n.) 


cogy/n 

A ’ 

q 2 gc 2 \'n+ 1 
A + J ’ 
sigVn + 1 
A 

q2S2gVn + 1 sigy/n+ 1 

A-k J A 
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Assuming all the frequencies/detunings are large com¬ 
pared to the Rabi frequencies of the drive fields, plugging 
these values into Hd and discarding all time-dependent 
terms, we have: 


H' = - 


Co/25 (j - 


aLcr -I- a'^aRC 


( 12 ) 


q2gc2f2 
' A + (5 


gfs 

2 


5252 

A-h<5 


a'j^a^e + a aRe^^ 
92S2 2si 


-b 


A + <5 


(«fl + afl) 


Aside from a spin-independent + ur^ term which 

can be canceled out by an oscillating voltage VQr sin ujRt 
applied to the resonator, this is precisely H' = 
A (h • ax) + clr^ as desired. This mechanism can 
thus be used to measure the x and y components of the 
qubits’ rotating frame states, allowing for dressed state 
measurement without having to turn off the many-qubit 
Hamiltonian. 

Note also that this formulation can be applied to cou¬ 
ple pairs of qubits, and generate the x, y, zx and zy com¬ 
ponents of CabO'iO' 2 - The zz components from 

this method, however, will be very small and are not 
easily tuned (they arise from the dispersive energy shift 
of each qubit based on the state of the other, and thus 
can only be tuned through adjusting the relative qubit 
detuning), so it is less versatile than the flux biased JJ 
of the previous section. A further complication arises 
when one wishes to engineer couplings of different types 
between a single qubit and many others; for the flux bi¬ 
ased JJ, the fact that the coupling is specified entirely by 
the flux through the loop between two qubits means that 
couplings between each of the two qubits and anything 
else are irrelevant to that term in the Hamiltonian. In 
contrast, for capacitively coupled qubits, the coupling is 
generated by driving an internal degree of freedom in the 
qubit itself, which generically will influence the couplings 
between that qubit and any other degrees of freedom it 
interacts with. Finally, we note that an additional qubit 
A can also be used for measurement through this mech¬ 
anism; if the pulse H' — A (h • ctt) is applied for a 
known interval then the operator h • ax can be measured 
simply by measuring 

Having outlined two routes to universal interactions 
between pairs of superconducting qubits, we will now dis¬ 
cuss two applications: passive cooling through frequency 
conversion, effectively replicating the NMR process of 
dynamic nuclear polarization in superconducting qubits, 
and higher-order interactions, where a single two-level 
qubit can be mapped to a n-level device (where the first 
n levels are linearly spaced and there is a large nonlin¬ 
earity to level n -b 1). Higher-order interactions can be 
used to engineer complex anyon models mniiMi with 
non-abelian exchange statistics, so they are of consid¬ 
erable interest to both condensed matter and quantum 
information physics. 


IV. DYNAMIC NUCLEAR POLARIZATION IN 
QUBIT ARRAYS 


As a simple application of the results of this work, we 
can use frequency-converting couplings as a dynamical 
cooling mechanism for groups of superconducting qubits. 
Specifically, we will demonstrate that by coupling a group 
of qubits to a second array of lossy qubits with signifi¬ 
cantly higher energies, we can engineer a passive cooling 
mechanism analogous to dynamic nuclear polarization in 
NMR [27], where magnetic polarization can be trans¬ 
ferred from electrons to nuclei, whose vastly smaller gyro- 
magnetic moment makes them otherwise difficult to po¬ 
larize at achievable temperatures. Our mechanism could 
also be used, with some restrictions, to cool many-body 
systems below ambient temperature (as was done in an 
earlier experiment using a protocol which was specific to 
flux qubits [28]), and thus could be useful for quantum 
simulation and adiabatic quantum computing. 

For simplicity, we consider identical and uncoupled 
primary qubits P, with system Hamiltonian H = 
o'|/2. In a thermal environment, the average spin 
(cr*) = — tanh (w/fesT), and we wish to cool the system 
to a lower effective temperature without actively measur¬ 
ing the system, increasing uj or decreasing the environ¬ 
ment temperature T. To do so, we introduce a second, 
“shadow” lattice of auxiliary qubits S Elm, which have 
excitation energies ujs and a fast relaxation rate Fs. To 
couple the two lattices, we introduce a frequency con¬ 
verting coupling as in Q, giving the total system Hamil¬ 
tonian 


H 



3 


+25 XI ^jS + ^3 cos At, 

3 


where A = wg — a; and we assume g <C w, A. As before, 
we work in a resonantly rotating frame by applying the 


unitary transformation !'!')—>• exp 


* (A/ 2 ) Y.j 


|vk). 


Our rotating frame Hamiltonian H' is then given by 


= X [-| (^3 + ^3s) + 9 + ^3 4s)}l4) 

3 


For the moment we consider g = 0. Modeling the 
environment as a thermally populated bath of harmonic 
oscillators with energies near w and us [29| , we note that 
at thermal equilibrium, the density of excited qubits is 
proportional to the ratio of the excitation and relaxation 
rates for the qubits to exchange excitations with the bath, 
so that 


P 


+ 



(15) 


Here, Fp = (1 -|- Nth) k and Fp = NthU, where k is the 
qubit-bath coupling and Nth = is the (appro¬ 

priately normalized) density of excitations at the relevant 
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energy in the thermal bath, which is exponentially small 
at low temperatures. 

On the other hand, the shadow qubit excitation den¬ 
sity pg will be proportional to ^-^sI^bT^ despite that 
the two groups of qubits are isoenergetic in (13). This 
is because, as described in [8] and elsewhere, when con¬ 
sidering continuously driven, open quantum systems, the 
bath degrees of freedom must be transformed in the same 
manner as the qubits. In this frame, both sets of bath 
modes now have frequency w, and the qubit-bath cou¬ 
plings + ^j(S)^(S)) unchanged, though 

any static couplings between qubits and bath modes at 
different energies will be rapidly oscillating and can thus 
be neglected. This is of course natural because moving to 
the rotating frame is just a unitary transformation used 
to simplify the analysis of the system’s dynamics, but the 
key point is that if A ^ fc^T, the shadow qubit popu¬ 
lation will be exponentially smaller and they can in turn 
act as a zero temperature bath for the primary qubits. A 
quick Fermi’s golden rule calculation of the two-step pro¬ 
cess where an excitation in the primary qubits tunnels to 
a shadow qubit and decays into the oscillator bath yields 
an additional contribution to the primary qubit loss rate 

r;: 


Tp —> Tp -|- 


4g^rg 
4g2 + r|- 


(16) 


As the equilibrium population is still given by the ratio 
of the excitation and relaxation rates, the new effective 
system temperature is given by 







4g^rg \ 
4g2+r|y'- 


(17) 


In realistic devices the induced decay rate can be two to 
three orders of magnitude larger than the system-bath 
coupling K, so while the reduction in effective tempera¬ 
ture is only logarithmic, it can still be substantial. One 
obvious use for this mechanism is to reduce the error rate 
in preparing ancilla qubits- most quantum algorithms re¬ 
quire many auxiliary qubits to be introduced at various 
stages in the algorithm, always in a known state, so a pas¬ 
sive cooling scheme could dramatically reduce the chance 
that an ancilla qubit is initialized in state |1) rather than 
|0) as intended, without having to perform the extra step 
of measuring the qubit directly. Of course, the passive 
cooling mechanism must be turned off while the qubit is 
being used for computation, but as the coupling is en¬ 
tirely generated by a drive field this is easy to do. 

Provided that the shadow qubit detunings can be ap¬ 
propriately controlled, this system can also be used to 
cool strongly interacting many-qubit arrays, with the 
shadow qubit energies tuned to match the energy cost of 
excitations about the ground state. However, we caution 
that for generic models and more general couplings than 
simple frequency conversion, the shadow qubits will in¬ 
duce a finite error rate of their own which is only quadrat- 
ically suppressed in the ratio of the couplings to the en¬ 
ergy cost, gjui. To see this, we consider the previous 


example, but replace the exchange coupling by a full xx 
term: 


H" 


^[~2 


(18) 


The g term now introduces an error channel, in which 
the XX spontaneously excites both qubits and the shadow 
qubit decays, leaving the primary qubit excited. Includ¬ 
ing this term, the effective temperature now reads: 


g ficj / fc T’gff 


jkuT I 


4g2r 

4g2+r; 


(19) 


For low enough temperatures (or large enough g/w) the 
induced error rate can be larger than the exponentially 
suppressed coupling to the outside world, and the shadow 
lattice can actually heat the primary qubits. One must 
therefore be careful when designing passive cooling sys¬ 
tems to appropriately suppress the induced error rate, ei¬ 
ther by ensuring that the coupling Hamiltonian can only 
correct errors (as in the simple case of independent spins) 
or by making the coupling weak enough that the induced 
error rate is minimal. Whether or not passive cooling 
below ambient temperature can be implemented success¬ 
fully thus depends on the detailed properties of the model 
and ratio of its excitation gap to the environment tem¬ 
perature. 


V. HIGHER-ORDER INTERACTIONS 
WITHOUT GADGETS 


As a second application of these couplings, we will now 
demonstrate how one can map superconducting devices 
to linear, n-level systems. For concreteness, we will con¬ 
sider the transmon qubit pair in FIG. 1(a); the two qubits 
have energy levels {0,0;^, 2wi — ai, 3uji — Pi...}. We fur¬ 
ther define: 


A = 012 - Wi; a^j = at - af, Ff- = \ni) (n^l. (20) 


For simplicity, we will consider the case where we wish 
to simulate bosons hopping between a pair of sites with 
local three-body interactions: 


H = coi + 0202^ — g (^a\a 2 + a\ai 

+ g (fli)^ + (®2) 


( 21 ) 


By default, there is a strong 2-body interaction (coming 
from the nonlinearity) which we must eliminate, and 
the qubits are detuned from each other so no tunneling 
occurs without an applied drive signal. To enable pho¬ 
ton tunneling and eliminate the 2-body interaction, we 
introduce the following drive signal for the flux bias F (t) 
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in ([^: 

F {t) = ^ + ko cos (At) (22) 

+ki cos ((A + ai) t) + k 2 cos ((A — a 2 ) t) 

+ks [cos ((A + ai 2 ) t) + cos ((A - 012) t)]; 


We now move to a rotating frame via the following trans¬ 
formation: 


14') —>■ exp 


i (^-aiPi + Aala2 - a2Pi'j t jvl/). (23) 


In this frame, the diagonal part of H is simply 
wi -I- 0202 ^ -I- UiPi + U2P2, where Ui = -Pf, 

the two-photon nonlinearities have been eliminated in 
the unitary transformation. If, as before, we identify 

sin pi = s (a| -|- Oi), the tunneling terms will be time- 


dependent, but for each tunneling term, there exists a 
component of \22\ which eliminates the time dependence. 
For example, the k^ cos At term removes the time de¬ 
pendence from jOI) (lOj -I- jlO) (Olj, ki cos ((A -|- ai) t) + 
cos ((A — 0:2) t) removes the time dependence from 
|02) (Ilj -H jll) (02| and |20) (llj -h jll) (20|, and so on. 
Provided that the energy spacings ai and A are all large 
compared to the net tunneling matrix element g, the 
rapidly oscillating terms can be ignored and we obtain 
the desired Hamiltonian (21). Further, we have achieved 
this without increasing the complexity of the physical 
circuit (as in perturbative gadget constructions such as 
[26]). 


As in all previous examples, the phases of the signals 
at different frequencies must remain locked to avoid gen¬ 
erating unwanted complex phases in the rotating frame 
tunnel couplings. Higher order interactions can be gen¬ 
erated by adding additional signal components to cancel 
out the higher nonlinearities (such as the 3-body nonlin¬ 
earity Pi). Likewise, the natural 2-body interaction can 
be reintroduced, with freely tunable magnitude and sign, 
by simply detuning the drive signals away from the ai. 
Note also that when counter-rotating terms and other 


higher order processes are taken into account, residual 
nonlinearities will likely remain, though for realistic pa¬ 
rameters they will be 2-3 orders of magnitude smaller 
than the bare nonlinearity of the device. 


VI. CONCLUSION 

We have demonstrated that a flux-biased Josephson 
junction can act as a universal interaction resource in 
the driven arrays of superconducting qubits, allowing ar¬ 
bitrary spin models to simulated dynamically by sim¬ 
ply tuning applied drive helds (and leaving the physical 
wiring unchanged). We have also shown how most of 
this flexibility can be replicated through simple capaci¬ 
tive couplings between qubits with internal flux lines, po¬ 
tentially simplifying the construction, and demonstrated 
that the frequency-converting nature of the couplings can 
be used to passively cool qubits below the temperature 
of their surroundings. Finally, we have shown how incor¬ 
porating multiple frequency components in a frequency¬ 
converting tunnel coupling can selectively eliminate qubit 
nonlinearities in the rotating frame, and thus gener¬ 
ate higher-order local interactions than simple 2-body 
terms. Our results thus considerably expand the “tool¬ 
box” for quantum simulation in superconducting qubit 
arrays, and demonstrate that extremely diverse classes 
of interacting spin models could be simulated with the 
same physical hardware. 


VII. ACKNOWLEDGEMENTS 

I would like to thank M. Hafezi, J. Koch, Y. Lu, V. 
Oganesyan, D. Schuster and S. H. Simon for useful dis¬ 
cussions. This material is based on work which was sup¬ 
ported by the Graduate Center of the City University 
of New York, EPSRC Grant Nos. EP/I032487/1 and 
EP/I031014/1, and Worcester College of the University 
of Oxford. 


[1] D. P. DiVincenzo, D. Bacon, J. Kempe, G. Burkard, and 
K. B. Whaley, Nature 408, 339 (2000). 

[2] E. Farhi, J. Goldstone, S. Gutmann, J. Lapan, A. Lund- 
gren, and D. Preda, Science 292, 472 (2001). 

[3] E. Farhi, J. Goldstone, S. Gutmann, and M. Sipser, 
arXiv:quant-ph/0001106 (2000). 

[4] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and 
S. Das Sarma, Rev. Mod. Phys. 80, 1083 (2008), 
URL http: //link, aps . org/doi/10.1103/RevModPhys . 
80.1083 

[5] A. G. Fowler, M. Mariantoni, J. M. Martinis, and A. N. 
Cleland, Phys. Rev. A 86, 032324 (2012). 

[6] K. C. Young, M. Sarovar, and R. Blume-Kohout, Phys. 
Rev. X 3, 041013 (2013). 


[7] M. Sarovar and K. G. Young, New J. Phys. 15, 125032 
(2013). 

[8] E. Kapit, J. T. Chalker, and S. H. Simon, 
arXiv: 1408.0959 (2014). 

[9] B. J. Brown, D. Loss, J. K. Pachos, G. N. Self, and J. R. 
Wootton, arXiv:1411.6643 (2014). 

[10] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008), URL http://link.aps.org/doi/10. 
1103/RevModPhys.80.885 

[11] R. J. Schoelkopf and S. M. Girvin, Nature 451, 664 
(2008). 

[12] M. H. Devoret and R. J. Schoelkopf, Science 339, 1169 
(2013). 

[13] M. H. Devoret, Quantum Fluctuations, Les Hourches, 






Session LXIII (1995). 

[14] M. R. Geller, E. Donate, Y. Chen, C. Neill, P. Roushan, 
and J. M. Martinis, arXiv: 1405.1915 (2014). 

[15] R. Barends, J. Kelly, A. Megrant, A. Veitia, D. Sank, 
E. Jeffrey, T. C. White, J. Mutus, A. G. Fowler, 

B. Campbell, et al., Nature 508, 500 (2014). 

[16] E. Zakka-Bajjani, F. Nguyen, M. Lee, L. R. Vale, R. W. 
Simmonds, and J. Aumentado, Nature Physics 7, 599 
( 2011 ). 

[17] A. Kitaev, Ann. Phys. 321, 2 (2006). 

[18] E. Kapit, Phys. Rev. A 87, 062336 (2013), URL http: 
//link.aps.org/doi/10.1103/PhysRevA.87.062336 

[19] M. Hafezi, P. Adhikari, and J. M. Taylor, Phys. Rev. B 
90, 060503(R) (2014). 

[20] J. Koch, A. A. Houck, K. L. Hur, and S. M. Girvin, Phys. 
Rev. A 82, 043811 (2010), URL http://link.aps.org/ 
doi/10.1103/PhysRevA.82.043811 

[21] E. Kapit, M. Hafezi, and S. H. Simon, Physical Review 
X 4, 031039 (2014). 

[22] J. M. Chow, J. M. Gambetta, A. W. Cross, S. T. Merkel, 

C. Rigetti, and M. Steffen, New J. Phys. 15 115012 
(2013). 


[23] Y. Chen, C. Neill, P. Roushan, N. Leung, M. Fang, 
R. Barends, J. Kelly, B. Campbell, Z. Chen, B. Chiaro, 
et ah, Phys. Rev. Lett. 113, 220502 (2014). 

[24] D. Solenov, S. E. Economou, and T. L. Reinecke, Phys. 
Rev. B 89, 155404 (2014). 

[25] J. Koch, T. M. Yu, J. Gambetta, A. A. Houck, D. 1. 
Schuster, J. Majer, A. Blais, M. H. Devoret, S. M. Girvin, 
and R. J. Schoelkopf, Physical Review A 76, 042319 
(2007). 

[26] E. Kapit and S. H. Simon, Phys. Rev. B p. 184409 (2013). 

[27] T. Maly, G. T. Debelouchina, V. S. Bajaj, K.-N. Hu, C.- 
G. Joo, M. L. Mak-Jurkauskas, J. R. Sirigir, P. C. A. 
van der Wei, J. Herzfeld, R. J. Temkin, et al., J. Chem. 
Phys. 128, 052211 (2008). 

[28] S. O. Valenzuela, W. D. Oliver, D. M. Berns, K. K. 
Berggren, L. S. Levitov, and T. P. Orlando, Science 314, 
1589 (2006). 

[29] C. Gardiner and P. Zoller, Quantum Noise: A Handbook 
of Markovian and Non-Markovian Quantum Stochastic 
Methods with Applications to Quantum Optics (Springer, 
2004) 


